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Abstract. We propose to use Stokes parameter as an entanglement witness 
for correlated EPR mixed states of light. Such states can be generated with 
a beam splitter acting on two mixed squeezed states of light. Stokes witness 
operators are closely related to the Hanbury-Brown and Twiss interference and 
can be used to test entanglement in balanced homodyne experiments involving 
fluctuations of quantum quadratures of the electric field. 
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1. Introduction 

In the famous paper questioning the reality of quantum mechanics, Einstein, Podol- 
sky and Rosen [ ^ introduced a quantum state with nonclassical correlations be- 
tween two spatially separated massive particles. In its original formulation this 
EPR correlated state was presented in the following general form 



^(xi,x 2 ) = 2J V'n (xi)u n (x 2 ). (1) 



In this remarkable formula we recognize the Schmidt decomposition of an entan- 
gled state. The EPR understanding of the probabilistic interpretation of quantum 
mechanics has been criticized by N. Bohr [ 2 . This famous debate about the founda- 
tions of quantum theory and the meaning of entanglement has drawn the attention 
of many physicists. 
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As an example of a correlated state given by Eq.JTJ, EPR have used the follow- 
ing wave function 

^{x 1 ,x 2 )~5{x 1 -x 2 ). (2) 

This wave function reflects the presence of perfect correlations between the coordi- 
nates of the two particles. If the first particle is at x\, the other one is at Xi with 
certainty. Similarly, there exists a perfect correlation in the momenta of the two 
particles 

*(p 2 ,p 2 ) ^ 5(pi +p 2 ). (3) 

The most important tools used in tests of quantum entanglement have been 
based on various Bell inequalities. These inequalities became a driving force behind 
a very active experimental branch of a quantum optics [ EH |Sj ■ Most of 
the experimental tests have been limited to entangled states of two qubits. For 
entangled systems described by continuous variables like the EPR wave function, 
Bell inequalities of some kind are violated as it has been demonstrated recently [|7|. 

The entangled states defined by the Schmidt decomposition are systems de- 
scribed by wave functions, i.e., are pure quantum states. For a quantum system 
described by a mixed statistical operator p, it turns out that the concept of quan- 
tum entanglement, as defined by QJ, has to be generalized. In the general case of 
a density operator, rather than a wave function, one uses the definition of quantum 
separability introduced by Werner A general quantum density operator of a 
two-party system is separable if it is a convex sum of product states, 

p = ^Pk p[ k) ® p% witn ^,Pk = 1 and p k > , (4) 

k k 

where p^ and p 2 k ^ are statistical operators of the two subsystems in question. If 
this criterion is not satisfied the state is called non-separable or simply entangled. 

The problem of establishing which mixed state is separable and which is not, is 
much more complex, and the question of experimental tests for such systems is still 
an open problem. Therefore, within the last decade a notion of an entanglement 
witness has generated a lot of interest as a experimentally feasible way of entangle- 
ment test for mixed states. The idea relics on constructing an operator which mean 
value is positive for separable and negative for entangled states 

. ., _ J > 0, p is separable, 
^ ' \ < 0, p is entangled. 

Such an operator has been discussed first for discrete variables [0. A general 
method for the experimental detection of entanglement for qubits is discussed in 
[ ll()j . The other experimentally feasible method relies on violating the entropy 
inequalities [ 111). Another entanglement witness for continuous variables, based 
on the measurement of the fluctuations of photocurrent difference between two 
correlated modes of detected light, has been introduced in [IT2"]. 
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A method relying on the Hanbury-Brown and Twiss interference, in which one 
can introduce an entanglement witness operator yy( ffBT ) has been presented in [ 

eg. 

In this paper we explicitly construct an operator which witnesses continuous 
variable entanglement of the EPR state. This operator can be expressed in terms of 
quantum uncertainties of the Stokes parameter. We show that the measurement of 
witness operator is equivalent to homodyne detection of quantum correlations and 
fluctuations of the electric field quadratures. If the value of the quantum uncertainty 
is below the shot noise the state is entangled. Our work focuses on general mixed 
gaussian states, involving as an example mixed EPR states. 



2. Thermal mixed squeezed states 

Let us start our discussion with a general one-mode mixed gaussian state. A single 
mode electric field operator, oscillating with frequency u>, can be written in the 
following form 

E = Eaiae-™* + aV"*) = y/2E Q {X x costut + X 2 smart), (5) 

where X\ = and X 2 = are amplitude and phase hermitian quadratures 

of the electric field. 

A one-mode gaussian state is completely characterized by the elements of its 
covariance matrix given by the following second-order moments 

(a'a) — n, (a 2 ) = —to, (6) 

where n is a mean number of photons in the mode, and to is a squeezing parameter, 
which we assume for simplicity to be real. We have assumed that (a) = (at) = 0, 
which can always be arranged with a suitable unitary shift of a and a' . 

Knowing those two parameters it is easy to construct a positive density operator 
for the state 

e -JV a f a-M"a 2 -Ma i2 
Pa = ^ ' (7) 

where N and M are known functions of n and m [ I14j . In this formula Z is the 
partition function. This mixed squeezed state leads to the following uncertainty of 
electric field quadratures 



AXi = \jn + i -to, AX 2 = ^Jn + i + to, (8) 
and the Heisenberg uncertainty relation is 

AX 1 AX 2 = J(n + ^) -to 2 , (9) 
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which imposes the following condition for the squeezing parameter 

< |m| < y/n(n + l). (10) 

It is clear that this condition guarantees that p a > 0. The notion that to is 
the squeezing parameter of quantum quadratures is obvious from these relations. 
Squeezing means in this case fluctuations below the vacuum noise, i.e., 

n < \m\. (11) 

For \m\ = y/ n(n + 1), this state reduces to the single mode pure state. 



3. Entangling photons 

One of the simplest way to entangle photons is to use a beam splitter. The action 
of the beam splitter is a nonlocal operation that may, in general, produce correlated 
photons. The use of a beam splitter to entangle squeezed light has been implemented 
in several experiments [ 1151 116| . This process is described as a 50/50 beam splitter 
(BS) transformation of a product state into a correlated state 

BS : p a ® p b ^ p cd . (12) 

The operators p a ,b are the density operators of the single mode mixed squeezed 
states with equal number of photons but opposite squeezing phases 

(a f a) = (tfb) = n, 

„2\ _ /j.2 



[rr 



(b 2 ) = -to. (13) 



The relation between input and output modes operators on a 50/50 beam split- 
ter have the well known form 

a+b a-b 

c = ^T> d =^T- (14) 

The two output beams form a correlated squeezed mixed EPR state with equal 
mean number of photons in each mode 

(etc) = (dU) = n (15) 

and with the following correlation parameter 

(a 2 ) - (6 2 ) 

(«V = g = ~ m ' ( ) 

It is clear that the beam splitter correlates the two beams 

BS : [n,m)a®{n,—m)b^(n,m)cd- (17) 
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For | to | = y/n(n + 1) this state reduces to the well known pure two- mode 
squeezed state, generated in a process of nondegenerate optical parametric amplifi- 
cation (NOPA) 



|NOPA) = £]V&»»>i 



(18) 



where p n = is a Bose- Einstein photon distribution. This state is the closest 

experimental realization, of the original EPR state. In the limit of infinite intensity, 
n — » oo, this state with p n ~ 1 becomes the EPR state © 



lim (xi,X2|NOPA) = y' (xi , X2\n, n) ~ 8(xi — x<£). 

71 — >00 — ' 



(19) 



n=0 



The separability condition for the two mode mixed state has been addressed 
in several publications. In this paper we follow the results obtained in [ 1171 I18j . 
Figure n shows the results computed for the general mixed EPR state. 
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Fig. 1. The values on and above the curve n = h (— 1 + vT-j- 4m 2 ) correspond 
to the physical states (with a positive density operator) - EPR mixed states. The 
states between the curve and the dashed line n = m are entangled. The EPR states 
above the dashed line are classically correlated, so separable. 



4. Witnessing entanglement 

In reference [ 113) it has been shown that the entanglement witness operator, related 
to the second-order Hanbury-Brown and Twiss interference, is given by the following 
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hermitian operator 



1 2a*<tfab + tf 2 a 2 + a^b 2 



W (HBT) = ±_ ™ » "»-r» u -ru u 

2 (: (J„ + 7 fc ) 2 :) 1 Uj 

The mean value of this operator allows to witness the entanglement in the outcoming 
modes c and d. Using Eq. (|14H we obtain 

Ti{W (HBT) p\ - ' m ' 2 - ^ ~ °' se P arable ' 



2(3n 2 + |to| 2 ) \ < 0, entangled. 
5. Hong-Ou-Mandel interference 

The entanglement witness which has been introduced in the last section is exper- 
imentally feasible quantity. It is directly related to the second-order visibility 
observed in the Hong-Ou-Mandel time resolved interference. 

The joint coincidence of detecting photons from mode c and d at delayed times 
t and t + t can be obtained directly from the following second-order temporal 
coherence function 

G (2) (t, t + r) = (f?H (i) (t + r) £7« (i + r) (t) ) . (21) 

Assuming the random phases in the output state and random phases typical 
for a stationary stochastic phase diffusion model in mode functions, we achieve the 
probability of the joint detection depending only on the delay [ I19j 



p(r) = l-v^ expj-(£) j, (22) 

where r c is a coherence time. For the EPR correlated state, the visibility of inter- 
ference fringes is equal to 

(2) _ n 2 + |m| 2 
V ~3n 2 + H 2 ' 
The visibility is directly related to entanglement witness 

(W (HBT) )= 1 _ v (2)_ (24) 

If the EPR state is separable, n > \m\, the visibility is v^ 2 ' < \ and the witness 
(I24|l is positive. For entangled states xft> > ~, the witness takes a negative value. 

6. Stokes Parameters 

The Stokes operators for the outgoing modes c and d are given by the following 
expressions 

cU + d^c etc - dU 
= ^ , b z = , (25) 
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S y = ^^,So = ^^. (26) 



These definitions of the Stokes operators are known as the Schwinger representa- 
tion of the angular momentum operators in terms of annihilation operators of two 
harmonic oscillators. These operators provide a very useful and efficient tool in 
the theoretical and experimental description of polarization entangled states [ I20j . 
In this paper we will use the Stokes parameters for continuous variable entangled 
squeezed states. 

The Stokes operators form a set of noncommuting operators 

[S x ,S v ]=iS z (27) 

and their joint measurement is limited by the Heisenberg uncertainty relation 

(AS v f(AS z f > |(^)| 2 , (28) 

which is equivalent to a positive density matrix condition n(n + 1) > |m| 2 . 

The entanglement witness (|20() can also be expressed in terms of the normally 
ordered Stokes parameters 

WHBTh = I ( : 5 " : ) (2 9) 

[ } 2 (:S 2 :) + (:5 2 :) + (:5 2 :) 1 > 



1 (jf) ~ gjjo) 



2 (S|) + (Sg) + (S2)-|(5 ) 



(30) 



The mean value of the witness takes a negative value if any of the following 
inequalities is obeyed 

(:S 2 X :) > (: S 2 :) + (: S 2 :), (31) 

(Si) > (S 2 ) + (S 2 Z )-±(S ). (32) 

As we shall see below, the measurement of the uncertainty of a single Stokes 
parameter can also constitute an entanglement witness. The measurement of AS Z 
seems to be the best choice for an experimental realization. This measurement can 
be done using balanced homodyne detection with a strong local oscillator beam a = 
\a\e llp . In Fig. [21 we have depicted an experimental scheme of such a measurement. 
The two incoming modes a and b Ijl3(l are correlated at a 50/50 beam splitter. 
The detectors measure intensities of the new displaced modes c and d. The Stokes 
parameter S z for the detected outgoing modes c and d is given by 

- (fid 

S z = C -±^Ll, (33) 

where the modes c and d are displaced by the strong pump 

c — c + a, d — d + a. (34) 
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Fig. 2. Balanced homodyne detection. The outgoing modes, c and d given by (|15|) 
and (|16|) . are combined with a strong local oscillator beam at a 50/50 beam splitter. 



Defining rotated quadratures as 

X C (<P) 



V2 



: 



the uncertainty of Stokes parameter S z is equal to 



(AS Z ) 2 



i(n+l) 



The above formula simplifies to 
(AS Z ) 2 = 



n(n + 1) — m 2 |a| 2 
2 + ~2 



^ L (l + 2(n + mcos2v?)) 



(35) 



(36) 



(37) 



where (X c ((p)Xd(^p)} = mcos2ip is the correlation between the two quadratures. 
Choosing tp = ^ and keeping quadratic terms in a only, the formula l|37|l reduces to 



(A5 z ) 2 ~^(l + 2(n-m)). 



(38) 



The above expression shows that if the EPR state is entangled, so if m > n the 
uncertainty of S z is below the shot noise 



(AS,) 2 < 



(39) 



In a similar way other Stokes parameters can also serve as an entanglement witness. 
Evaluating the variance of the Stokes parameter S x 



S x = 



tfd + cftc 



with (S x ) = |a| 2 we get 
(AS*) 



(40) 



(41) 
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In the strong field limit this formula simplifies to 

(AS X ) 2 ~ i^- (l + 2(n-mcos2<p)). (42) 

Choosing (p = and keeping quadratic terms in a only, the formula (|42|l reduces to 
formula (|38[> . The outcoming state is entangled if (A5 X ) 2 is smaller than for a shot 
noise. Similarly for S y 

Its mean value (S y ) = and the variance is equal to 



Sy = 1^±. (43) 



^S y )^ n ±±4±^ + ^({xA V + ^-XA^^y). (44) 



This formula reduces to formula l|42|) giving the same condition for non-separability 
of outcoming EPR mixed state. 



7. Conclusions 

We have proposed a measurement of the uncertainties of Stokes parameter S = 
(S x , S y ,S z ) as an entanglement witness for a general class of EPR mixed correlated 
states. This method seems to be relatively easy to implement experimentally If the 
quantum fluctuations of the Stokes paramaters is below the shot noise (the value 
obtained for a vacuum state), the state is entangled i.e., non-separable. 
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